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ABSTRACT 


The  size  of  the  shock-layer  governed  by  a  conservation  law  is  studied.  The  conservation 
law  is  a  parabolic  reaction-convection-diffusion  equation  with  a  small  parameter  multiplying 
the  diffusion  term  and  convex  flux.  Rigorous  upper  and  lower  bounding  functions  for  the 
solution  of  the  conservation  law  are  established  based  on  maximum-principle  arguments. 
The  bounding  functions  demonstrate  that  the  size  of  the  shock-layer  is  proportional  to  the 
parameter  multiplying  the  diffusion  term 


*  Research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract 
No.  NAS1-1S605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 


V 


1.  Introduction.  Upper  and  lower  bounding  functions  are  presented.  The  bounds 
demonstrate  that  solutions  to  the  singularly  perturbed  hyperbolic  partial  differential  equa¬ 
tion 


nf  ,  du  df(u)  d2u 
(1)  + 

with  a  pre-existing  shock  have  shock-layers  of  width  0(e).  The  analysis  is  performed  in  the 
style  of  Howes  [4,  5,  6).  It  begins  with  a  multiple-scales  asymptotic  analysis.  This  provides 
the  appropriate  local  scalings  and  indicates  candidate  forms  for  a  bounding  function.  A 
bounding  function  is  constructed  from  these  candidates.  Maximum  principle  arguments  are 
then  used  to  rigorously  establish  bounds  tor  the  solution.  In  this  way,  upper  and  lower 
bounds  for  solutions  to  equation  (1)  are  established  for  any  e. 

This  presentation  will  concentrate  on  a  comparison  between  solutions  to  (1)  and  solutions 
to  the  corresponding  reduced  equation 


(2) 


+ 


d_m 

dx 


-  R(U)  =  0 


obtained  by  setting  e  =  0.  The  result  will  be  a  bound  on  the  difference  between  the  solution 
to  this  reduced  equation  and  the  solution  to  equation  (1). 

There  are  implications  of  this  analysis  for  the  computational  aspects  of  the  problem 
as  well  as  for  the  physics  modeled  by  conservation  laws.  The  bounding  functions  result  in 
an  upper  bound  on  the  size  of  the  shock-layer.  They  isolate  the  internal-layer  region  in 
which  viscosity  is  important  from  the  convection-dominated  outer  region.  This  validates  the 
assumption  that  the  lack  of  resolution  of  the  physics  in  the  shock-layer  effects  the  solution 
in  smooth  regions  very  little  under  certain  circumstances;  thus,  the  hyperbolic  equation  (2) 
may  be  substituted  for  equation  (1). 

This  work  sharpens  the  bounds  of  Howes  [5,  7]  which  apply  in  a  more  general  setting. 
The  shock-birth  region  will  not  be  studied  here. 


2.  Problem  Specification.  This  paper  is  concerned  with  presenting  bounds  for  the 
solution  u  to  the  quasilinear  parabolic  equation  (1)  on  the  domain 


(3) 


D  :=  {(s,t)|0  <  x  <  b,0  <  t  <  T}, 


subject  to 


(4)  u(x,0)=^(x),  0  <  x  <  6; 

(5)  u(0,t)  =  a(t),  Q<t<T\  and 
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(c)  u(b,t)  =  P(t),  0  <t<T. 

Let  the  portion  of  the  boundary  along  which  the  data  is  specified  be  denoted  by 


n  :=  {(x,t)|0  <x<b,t  =  0}(J{OM)|0  <  t  <  T,  x  =  0,6}. 

The  boundary  and  initial  data  are  continuous  and  sufficiently  smooth  so  that  u  is 
uniquely  defined  [1].  In  particular,  the  boundary  data  a  and  f3  have  derivatives  bounded 
independent  of  e.  Corner-layers  in  u  are  prohibited  by  assuming  the  compatibility  conditions 

(7)  a(0)  =  5(0),  g(b)  =  /?( 0) 

(8)  ^  +  ^f(g)  -  R{g )  =  o,  for  (x,  t)  =  (0,0); 
and 

(9)  ^  +  ^/(<0  ~  =  °>  for  *)  =  °)' 

For  simplicity,  it  is  assumed  that  all  boundaries  are  inflow  boundaries;  hence,  ct(t)f'(a(t))  > 
0  and  <  0.  Also  assume  that  there  is  a  single  shock-layer  in  the  initial  data  that 

is  contained  in  an  0(eln(e))  neighborhood  tt0  of  (Fo,  0),  where  To  is  the  location  of  the  shock 
in  U  at  t  =  0.  The  domain  of  the  initial  viscous-layer  is  7r0  :=  {(x,t)|t  =  0  and  |x  —  r0i  < 
dn(e)}. 

The  solutions  to  the  parabolic  problem  will  be  compared  to  a  weak  solution  of  the 
hyperbolic  equation  (2).  Let  U  be  the  weak  solution  of  (2)  with  boundary  data  (5-6)  that  is 
the  solution  to  (1)  in  the  limit  as  e  >  0  tends  to  zero  (denoted  as  e  |  0).  The  initial  condition 
will  reflect  a  shock  eminating  from  (t,x)  =  (0,ro).  Thus,  the  initial  condition  for  U  is 

(10)  U(x,0)  =  g0(x),  0  <  x  <  b, 

where  the  difference  g  —  50  is  zero  except  in  7T0.  The  relationship  between  5  and  g0  in  7r0  will 
be  discussed  in  more  detail  in  the  proof  of  Theorem  4.4. 

Let  the  path  of  the  shock  in  U  be  given  by  the  curve  (x,<)  =  (F(t),f).  It  is  natural  to 
describe  the  values  of  U  at  the  shock  as 

Ur( t)  =  lim  U(x,t) 
xir(f) 

and 


17l(0  =  lim  U(x,t). 
itr(f)  v  ' 
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For  ease  of  presentation  we  will  assume  that  Ui>  UR.  The  solution  U  to  (2)  will  satisfy  the 
entropy  condition 

(11)  UL(t)  >  5(f)  >  UR{t) 

where  the  speed  S(t)  of  the  shock  is  given  by  the  Rankine-Hugoniot  jump  condition  [8] 

_  mm  -  mm 

1  }  (  )  uL{t)  -  uR(t)  • 

The  entropy  condition  may  be  written  as 

(13)  n(t)  =  UL{t)  -  UR(t)  >  fio, 

for  a  constant  >  0  that  is  independent  of  e.  It  is  assumed  that  the  shock  is  part  of  the 
initial  data  and  exists  for  the  entire  domain  considered;  thus,  /i(f)  is  defined  for  all  t  >  0. 
Notice  that  the  initial  condition  for  T  may  be  obtained  from  g0  and  that  F  satisfies  the 
ordinary  differential  equation  dV(t)/dt  =  S(t).  Thus,  T  is  uniquely  defined. 

3.  Asymptotic  Analysis.  Some  of  the  relevant  physics  is  presented  here.  The  discus¬ 
sion  includes  introducing  the  appropriate  scales  and  equations  when  equation  (1)  is  used  to 
model  a  shock-layer.  These  results  are  obtained  using  a  heuristic  analysis  that  treats  e  as  a 
small  parameter.  The  results  are  exploited  in  the  construction  of  the  bounding  function  and 
are  made  more  rigorous  in  Section  4.  First  some  of  the  properties  of  the  inviscid  solution 
will  be  discussed. 


3.1.  Derivation  of  Bounding  Function.  In  this  section  we  will  derive  a  canonical 
form  for  the  bounding  function.  The  assumptions  utilized  in  the  derivation  are  used  only  to 
motivate  the  equation  governing  the  canonical  form.  The  assumptions  are  not  necessarily 
valid.  The  validity  of  the  bounding  function  will  be  established  when  it  is  used  in  Section  4. 

The  canonical  form  for  the  solution  of  the  viscous  problem  is  derived  here.  Large  gra¬ 
dients  in  the  neighborhood  of  a  shock-layer  are  resolved  by  using  the  spatial  scale  1/e.  This 
scale  is  combined  with  shock-following  to  obtain  the  internal-layer  coordinates 


x  -  r 

(  =  -  an(J  T  =  t. 


The  transformation  defined  by  these  coordinates  is  applied  to  equation  (1)  to  obtain 


d/(£)  _  c(A^.  _  ^ 

d(  [>d(  dC2 


where  =  u(x,t).  This  suggests  the  regular  expansion 


u  =  fta  +  eux  +  .... 
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for  it.  We  assume  that  this  expansion  is  a  •priori  valid  in  the  shock  layer  and  use  identification 
in  e  to  obtain  the  equation 


dfjuo) 

d( 


v duo 

d( 


d2ito 

a? 


for  Uq.  This  equation  is  integrated  with  respect  to  (  once  to  obtain 

/(uo)  -  S(t)iio  -  =  const. 


To  make  this  equation  easier  to  solve,  we  approximate  the  term  S  —  f/ito  by  k(p,  and  we 
set  the  constant  in  the  right-hand-side  of  the  above  equation  to  zero.  The  solution  to  this 
equation  is  the  exponential 


(14) 


exp(— k(p+1). 


Since  the  magnitude  of  the  bounding  function  should  be  decaying  away  from  the  shock- 
layer,  we  assume  k  >  0,  and  substitute  |£|  for  (.  This  is  the  candidate  form  for  the  bounding 
function. 

Guided  by  the  form  (14),  let  the  general  form  for  the  bounding  function  u  be 
(15)  ui  =  ai/t  -f  6e~KiP . 


Here  i:(r(t),  x,  e)  is  a  linear  measure  of  the  distance  between  (a :,t)  and  {T(t),  t).  The  power 
of  the  distance  function  is  some  positive  p,  and  0  is  a  non-negative  function  of  t. 


4.  Comparison  Theorem.  In  this  section  we  will  show  that  the  function  w  +  U  is  an 
upper  bound  for  u.  The  statement  of  this  result  is  in  the  form  of  a  comparison  theorem. 
Several  lemmas  needed  in  the  proof  of  Theorem  4.4  will  be  presented  first.  The  first  lemma  is 
a  maximum  principle  stated  in  the  form  most  useful  for  the  proof  of  the  theorem.  Lemma  4.2 
demonstrates  how  to  choose  the  parameters  in  w  so  that  the  maximum  principle  is  satisfied 
within  each  of  the  regions  Q0{t)  '■=  {®|0  <  x  <  F(t)}  and  Qfft)  :=  {m |F(t)  <  x}.  Finally,  we 
see  how  the  parameters  in  w  can  be  chosen  to  satisfy  the  maximum  principle  on  (x,  t)  —  (T,  t). 

4.1.  Maximum  Principle.  This  lemma  is  a  modification  of  the  the  Nagumo-Westphal 
Lemma  [11]  to  include  functions  that  are  C00  except  on  a  set  of  measure  zero  where  they 
may  be  only  C°.  A  condition  on  the  spatial  derivative  replaces  the  condition  involving  a 
parabolic  operator  on  this  set.  This  lemma  is  a  direct  extension  of  the  result  by  Nagumo 
and  Westfal;  thus,  it  is  presented  without  proof. 


Lemma  4.1.  Let  z(x,t)  be  a  continuous  function  that  is  differentiable  except  on  a  finite 
number  of  curves.  Suppose 


(16) 


P[z}  >  0 
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in  the  regions  where  z  is  differentiable,  while  z  satisfies 

(i7)  >  ’Js^r  •  mw 


for  curves  ( x,t )  =  (X(t),t)  on  which  z  is  continuous  but  not  differentiable.  When  these 
conditions  are  satisfied  along  with 


z  >  u  for  (x,t)  €  II, 


then  z>u  throughout  D. 


The  implications  of  this  lemma  at  discontinuities  on  the  choice  of  the  bounding  function 
are  demonstrated  in  Figure  1.  The  jump  in  the  first  derivative  with  respect  to  x  is  larger  on 


Fig.  1.  Discontinuous  Upper  Bounding  Function 
the  left  of  ( x,t )  =  (X,t)  than  on  the  right. 

4.2.  Analysis  of  co:  Continuous  Regions.  The  results  in  this  section  require  a  con¬ 
vex  flux  function.  Let  Q(v)  =  The  assumption  of  convex  flux  means 

Q'(v)  >  const 

for  all  v  and  for  a  positive  constant  independent  of  e  where  Q'(v)  —  We  also  assume 
that  Q'{y )  is  bounded  above. 

The  shock  speed  S{t)  and  the  values  of  Q(u)  are  related.  The  mean  value  theorem  states 

that 

(19)  S(t)  =  Q(U) 
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for  some  intermediate  function  U(t)  in  the  open  interval  \Ur,Ul[.  We  will  assume  the 
slightly  more  restrictive  case  of  U  in  the  closed  interval  [Ur  +  A,  Ui  —  A],  where  A  is  a 
positive  constant.  This  relationship  between  S  and  Q  will  be  exploited  in  the  proof  of  the 
following  lemma. 


Lemma  4.2.  There  is  an  close  enough  to  F  and.  an  &q  small  enough,  such  that  for 
each  region  Q,  we  may  construct  a  specific  form  of  uj  to  satisfy  inequality  ( 16 ). 

Proof.  First  some  algebraic  details  are  discussed.  Let  x  =  \x  —  Y\/v.  Differentiation  of 
u  results  in 


=  av 


6' 


dpKS 


xp~ 


-Ki? 


and 


wx 


QpK  -p_i  -Kir 

- x y  e  , 

v 


+  c.K„ 


We  use  the  mean  value  theorem  to  obtain  f(U  +  u>)x  =  Q(U  +  uj)ujx  +  (Q(U)+ioQ')Ux,  where 
Q'  is  evaluated  at  some  function  between  U  and  U  +w.  The  parabolic  operator  applied  to 
the  bounding  function  is 


P[U  +  w]  =  wt  +  Qux  +  Q'uUx  -  e.Uxx  -  -  R(U  +  w) 


or 

P[U  +  w]  =  er!  +  r2  +  —rTs  +  -U, 
v2  v 

where 


rx  =  a  -  Uxx  +  atQ'Ux , 
t2  =  [O'  +  8QUx)e~Kir  -  R, 
t3  =  -8K2p2x2^  ]e'KiP 


and 


t4  =  8Kp[{S  -  Q)xp~l  +  (p  - 
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We  will  consider  two  cases  based  on  the  two  forms  of  uj  to  verify  that  (16)  is  satisfied: 
Case  I.  When  x  £  Q0  we  lot  u  =  uq,  where 


ttio  =  o.et. 

We  choose  a  large  enough  and  T  small  enough  so  that  Tj  >  0.  Inequality  (16)  holds  in  this 
case  since  t2,  T3,  r4  =  0. 

Case  II.  When  i  £  fli  we  let  uj  =  wlt  where 

u>i  =  aet- f  . 

Here,  6  and  k  are  positive  constants  independent  of  e,  and  p.  is  a  positive  function  of  t.  The 
term  Ti  is  positive  from  the  choices  in  Case  I;  thus,  it  is  sufficient  to  consider  only  r2,  r3  and 
m-  We  choose  v  >>  e.  For  example,  we  could  use  v  —  e\  ln(e)|.  Then,  we  assume  that  e0  is 
small  enough  that  With  this  choice  of  1/,  it  is  sufficient  to  show  that  r4  is  bounded  below  by 
a  positive  constant  that  is  independent  of  e. 

First  we  will  establish  a  lower  bound  for  Sit)  —  Q(U  +  oj).  These  properties  are  based 
on  the  fluid  dynamics  properties  of  the  problem.  Namely,  Q(u )  is  the  speed  at  which  the 
characteristics  of  u  travel.  The  characteristics  will  be  traveling  faster  than  the  shock  for 
x  to  the  left  of  T,  and  slower  than  the  shock  for  x  to  the  right  of  F.  This  means  that 
5(f)  —  Q(U  +  w)  is  positive  for  x  >  &o(/c).  The  relation  between  x  and  S  —  Q  is  depicted 
in  Fig  2.  (The  precise  shape  of  the  curve  cannot  be  known  without  more  information  about 
Q.)  The  location  of  b0(n)  moves  closer  to  the  origin  as  k  increases.  Also,  multiplying  S  —  Q 
by  x  reduces  the  variation  near  the  origin.  It  is  clear  that  we  can  choose  k  large  enough  such 
that 

1  +  (5  -  Q)x  >  i. 

□ 


We  now  have  specific  forms  for  the  bounding  function  in  each  of  the  domains.  Next  we 
must  show  that  we  can  satisfy  inequality  (17)  and  C°  continuity  with  U  +  w  with  the  scaling 
and  k  of  the  previous  lemma. 

4.3.  Analysis  of  w:  Discontinuities.  First  we  will  establish  the  C°  continuity  by 
choosing  the  parameters  for  the  bounding  function.  Condition  (18)  will  also  be  satisfied  by 
the  choices  in  this  lemma. 

Lemma  4.3.  We  may  choose  the  parameters  in  Wj  to  satisfy  (18)  and  such  that  w  is 

C°. 
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Fig.  2.  Example  form  for  S-Q 

Proof.  To  obtain  C°  continuity  let  Q  —  p. 

The  result  will  follow  providing  (18)  is  satisfied  for  (x,t)  £  II.  This  relation  clearly  holds 
except  possibly  for  (x,f)  £  ~a-  In  this  region  it  is  necessary  to  impose  some  restriction  on 
g  —  go .  Assume  that  the  viscous  profile  of  the  shock-layer  is  of  an  exponential  type  that  is 
bounded  by  the  function  ui  presented  in  this  proof.  This  is  not  a  severe  restriction,  since  a 
layer  with  a  profile  like  that  of  tanhf— (cc/e)]  (a  solution  to  Burgers’  equation)  is  suitable.  □ 

Theorem  4.4.  Assume  that  g  and  go  with  their  first  and  second  derivatives  are  bounded 
independent  of  e  except  that  g0  has  a  jump  at  (x,t)  =  (ro,0)  and  g  may  have  derivatives 
bounded  depending  on  e  for  (x,t)  £  tt0.  Suppose  w  is  constructed  fromui0  and  uq  of  Lemma 
4.2.  Also  assume  that  Ux  is  continuous  across  (a;,  t)  =  (T,  t).  Then  there  is  a  positive  eo  such 
that 

(20)  U  +  w  >  u 

for  ( x ,  t)  £  D  when  eo  >  e  >  0. 

Proof.  We  will  have  only  one  case  to  consider  based  on  imposing  (17)  on  the  curve 
(x,t)  —  (r(f),f).  Consider  a  modified  distance  function  x  =  |*  —  T  +  6v\/v.  Observe  that 
for  positive  5,  the  term  wx  from  equation  (4.2)  is  negative  and  has  magnitude  0(5).  If  we 
take  6  asymptotically  close  to  zero,  then  we  still  satisfy  this  condition.  Thus,  by  taking  the 
limit,  we  may  return  to  the  original  basis  function  since  lim^  i0  x  =  x.  □ 

Remark  2.  The  constraint  on  Ux  across  (x,t)  —  (T,t)  can  be  relaxed  somewhat.  This 
constraint  was  imposed  so  that  S  could  be  taken  asymptotically  close  to  zero  in  the  proof. 
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With  a  more  careful  choice  of  8  this  constraint  can  be  eased. 


Corollary  4.5.  Under  the  conditions  of  the  theorem, 

k\x  —  r(t)| 


(21)  —  I/|  =  0  ^/iexp 

when  £q  >  e  >  0  and  v(e)  is  chosen  as  above. 


+  0(e), 


This  result  follows  directly  using  symmetric  arguments  to  obtain  a  lower  bound. 


5.  Implications.  A  direct  result  of  this  theorem  is  an  upper  bound  on  the  size  of  the 
shock-layer.  In  this  context,  the  shock-layer  is  defined  as  the  region  in  which  the  solution  to 
(1)  differs  from  the  solution  to  (2)  by  more  than  a  specified  amount.  Namely,  it  is  the  region 
in  which 


(22)  \u-U\>o 

for  some  positive  a.  As  reflected  in  the  following  corollaries,  there  are  different  results 
depending  on  whether  a  is  a  function  of  e  or  not.  This  result  follows  directly  from  the 
theorem.  The  following  corollary  extends  this  result  to  the  case  when  c  is  independent  of  e 
and  is  a  direct  result  of  Corollary  4.5. 

Corollary  5.1.  Suppose  the  conditions  of  Theorem  4.5  obtain.  Let  (x(t),t)  be  the 
independent  variables  for  which  Inequality  (22)  is  satisfied.  If  a  is  a  constant  independent  of 
e  then  there  is  an  e0  small  enough  so  that  \u  —  U\  <  o  when 

\x(t)-m\<o(i') 

for  e0  >  e  >  0. 

When  we  define  the  shock-layer  as  the  region  such  that 

|u(x,t)  —  U(x,t)\  >  o, 
then  there  is  a  positive  eo  such  that 

|x(f)  -  T(<)|  <  0(eln(e)) 

for  e0  >  e  >  0. 

Physically  motivated  domain  decomposition  algorithms  can  be  based  on  ideas  presented 
in  the  analysis  discussed  herein.  The  computational  domain  can  be  partitioned  into  subdo¬ 
mains  that  have  different  physical  behavior.  Inside  each  of  these  regions  different  modeling 
equations  (and  hence  different  numerical  methods)  are  used  so  computational  effort  is  con¬ 
centrated  on  the  relevent  physics  [10],  [2],  Global  error  bounds  have  been  developed  for  these 
domain  decomposition  methods  [9],  In  addition,  the  general  idea  of  using  different  modeling 
equations  can  be  extended  to  systems  of  equations  [3]. 
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